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We report on the calculation of the cyclotron spin-flip excitation (CSFE) in a 
spin-polarized quantum Hall system at unit filling. This mode has a double-exciton 
component which contributes to the CSFE correlation energy but can not be found 
by means of a mean field (MF) approach. The result is compared with available 
experimental data. 
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A two-dimensional electron gas (2DEG) in a high perpendicular magnetic field possesses 
many remarkable features.- In particular, it presents a rare case of strongly correlated system 
governed by real Coulomb interaction (not by a model Hamiltonian!) where, nevertheless, 
some solutions of the quantum many-body problem can be found exactly. Indeed, under the 
conditions of integer quantum Hall effect (when the filling factor is v = 1, 2, 3, ...), the one- 
cyclotron magnetoplasma and the lowest spin-flip modes are calculated analytically to the 
leading order in the parameter r c = Ec/tiLV c .^^ [u c is the cyclotron frequency; Eq = ae 2 / kIb 
is the characteristic interaction energy, a being the average form-factor related to the finite 
thickness of the 2DEG (0.3 < a < 1); Is is the magnetic length.] This astounding property 
is the feature of either filled or half-filled highest-occupied Landau level (LL) where the 
simplest-type excitations are single excitons or superposition of single-exciton modes. The 
many-body problem is thereby reduced to the two-body one, i.e. to the interaction of 
electron with an effective hole. Being quite in the context of similar studies, the present 
letter concerns however the case which can not be reduced to a single-exciton problem. 

We remind that 2DEG excitons are characterized by sublevels a — (n a , o~ a ) and b = (n&, cr&) 
where electron is promoted from the n„-th LL with spin-component S z = a a to the n^-th 
LL with S z = ab. The relevant quantum numbers are 5n = rib — n a , SS Z = ab — cr a) and 
the two-dimensional (2D) wave vector q. The single exciton problem is exactly solvable 
in the following cases: (i) at odd filling v when 5n = 1 and 5S Z = (magnetoplasmon) 
or Sn = and SS Z = —1 (spin wave);-^ (ii) at even v when 8n = 1 and d~S z = 0, ±1 
(magnetoplasmon and spin- flip triplet).- 1 ^ At the same time the two-body problem may be 
discussed within a MF approach (in some publications called 'time-dependent Hartree-Fock' 
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approximation-^) which excludes any quantum fluctuations from a single exciton to double- 
or many-exciton states. For the above simplest cases of Sn and SS Z , the MF calculation gives 
an asymptotically exact result which may be found perturbatively to the first order in r c ,- 
because these (Sn, SS Z ) sets can not correspond to any states except single-exciton modes. 
Any complication of (Sn, SS Z ) makes the calculations substantially more difficult due to the 
necessary expansion of the basis to the entire continuous set of many-exciton states with the 
same total numbers Sn, SS Z , and q. For example, the double-cyclotron plasmon with Sn = 2, 
SS z = and with given q 'dissociates' into double-exciton states consisting of one-cyclotron 
plasmon's pairs with the total momentum equal to q. 4 At odd v, a similar 'dissociation' 
occurs for the CSFE, where Sn = —SS Z = 1. The proper double-exciton states are pairs of 
a magnetoplasmon (Sn = 1, SS Z = 0) and a spin wave (Sn = 0, SS Z = — 1). The problem 
thus changes from the two-body case to the four-body one, and the correct solution should 
be presented in the form of combination of the single-exciton mode and continuous set of 
double-exciton states. 6 It is important that in both cases the desired solution corresponds 
to a discrete line against the background of a continuous spectrum of free exciton pairs. 
The technique of correct solution has to be of essentially non-Hartree-Fock (non-HF) type. 
Actually this letter concerns the fundamental question of consistency of the MF approach. 

By considering the case of unit filling factor where the number of electrons is equal to the 
number of magnetic flux quanta N^, now we report on a study of the CSFE with q = 0. This 
state is optically active and identified in the ILS experiments.— 1 ^ 1 - Besides, it is exactly this 
spin-flip magnetoplasma mode which is the key component of the elementary perturbation 
used in the microscopic approach to the skyrmionic problem.— The calculation is performed 
in 'quasi-analytical' way which should, in principle, lead to the result which is exact in the 
leading approximation in r c . In our case the envelope function determining the combination 
of the double-exciton states is one-dimensional — i.e., it only depends on the modulus of the 
excitons' relative momentum. This function is chosen in the form of expansion over infinite 
orthogonal basis, where every basis vector obeys a specific symmetry condition necessary for 
the total envelope function. Even to the first- order approximation in r c , we obtain a non-HF 
correction to the former HF result-* 1 ^ for the CSFE energy. 

As a technique, we use the excitonic representation (ER) which is a convenient tool for 
description of the 2DEG in a perpendicular magnetic field.-^ 1 ^ When acting on the vacuum 

|0) (in our case |0) = | f , t, ••• t )), the exciton operators produce a set of basis states which 
diagonalize the single-particle term of the Hamiltonian and some part Hed of the interaction 
Hamiltonian.-* 1 ^ Exciton states are classified by q, and it is essential that in this basis the 
LL degeneracy is lifted. 
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So, the generic Hamiltonian is H = Hi + H mt where 



^ (iV -eA/ c y+gfi B BS z 



^n-fr) and 



(1) 



01,(72 



Choosing, e.g, the Landau gauge and substituting for the Schrodinger operator W a = 
^np^npa^npa (indexes n,p,<T label the LL number, intra-LL state, and spin sublevel), one 
can express the Hamiltonian (1) in terms of combinations of various components of the 
density-matrix operators.-'&i 2 - These are exciton operators defined as—&^^ 



^abq 



N <t>~ 1 ' 2 £ e ~ iqXPb 



P+ 



Qy_ a p-^ 



and Q afeq = Q\ a _ c 



(2) 



and obeying the commutation algebra ( 



-*^cd qi ' ~^ab q2 



= N 



-1/2 



( „-i(qiXqa)«/2 X q] _ i(qiXq 2 ) z /2 r £)t 

fc,c^ ac ; qi+q2 c °a,d^cb qi+q 2 
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(in our units /# = W ch/eB = 1). Here a, 6, c, ... are binary indexes (see above), which means 



= % ... We will also employ for binary indexes the notations n= (n, f ) 



that a p — a nap(Ja . 

and n = (n, |), so that the single-mode component of the CSFE is defined as Q} Q j |0). The 
interaction Hamiltonian can be presented as H mt = Hed + H' where Hed , if applied to the 
state Ql bq \0), yields a combination of single-exciton states with the same numbers Sn, 5S Z , 
and q (see Refs. |f3 12 and therein Hed expressed in terms of exciton operators). In the 
framework of the above HF approximation, the CSFE correlation energy 8,10 is obtained 
from the equation £ j(q) = (0| QoTq[-^int, Q^j ]|0) where only the Hed part of the interaction 
Hamiltonian contributes to the expectation. In the following, we need this so-called HF 
value at q — 0, namely £ j(0) = £hf = § p 3 dpV(p)e~ p2 ^ 2 where 27rV(q) is the Fourier 
component of the effective Coulomb vertex in the layer.— (In the strictly 2D limit a — > 1, 
and V(q) — > e 2 / ' nlsq) 

The problem arises due to the 'troublesome' part H' of the interaction Hamiltonian which 
can not be diagonalized in terms of single-exciton states. For our task we keep in H' only the 



terms contributing to 



H' Q, 



Olq 



|0) and besides preserving the cyclotron part of the total 



energy (i.e. commuting with Hi). In terms of the ER these are 



5S 



^Ev^'^oiqQoTq + H.c. 



(4) 



Using Eqs. (3) and identities QtJO) =N~ l/ X,o|0> if a=(0,T) and QLjO)=0 if a^(0,T), 
one can find that the operation of H' Q j on vector Qqj |0) results in a combination of states of 
the type of A^ _1//2 £ s f(s)Q^ / 2 _ s Qo lq / 2 4s|0) with a certain regular and square integrable 
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envelope function, J|/(s)| 2 <is~ 1. The norm of this combination is not small as compared 
to (0| 2oTq2oi q |0) = 1) an d the terms (4) must be taken into account when calculating the 
CSFE energy. 

On the other hand, if the set of double-exciton states |s, q) = Qoo q / 2 -s2oiq/24slO) 
is considered, then one finds that they, first, are not exactly but 'almost' orthogonal: 
(qi,si|s 2 ,q 2 ) = <5 qi , q2 {<5 Sl ,s 2 }, where {5 Sl , S2 } = 5 SuS2 - e l(siXS2)z /^; and, second, |s, q) sat- 
isfies the equation 



flint, Q 



00q/2-s^01q/24s 



(5) 



|0> = [£ 8w (|q/2 + s|) + £ mp (|q/2-s|)] |s,q) + |e>, 

where the state \e) has a negligibly small norm: (e\e) ~ Eq/N^. Therefore the double- 
exciton state |s,q) in the thermodynamic limit actually corresponds to free noninteracting 
excitons: one of them is a spin exciton (spin wave) with energy \gfisB\ +S SW where 

pd P V(p)e~ p2 / 2 [1-J (pq)], 



£ S w(q) 



while the other is a magnetoplasmon with energy hu c +S mp where 



-9 2 /2 



poo 

+ / pdpe' 



p2/2 v( P ) 1 



p- 
~2 



(6) 
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[Jo is the Bessel function (cf. Refs. 

Thus we try for the CSFE state the vector |X q ) = X q |0) where X q is a combined operator 



(8) 



Actually only a certain 'antisymmetrized' part {^p q } of the envelope functions contributes 
to the double-exciton combination in |X q ).-^ In our case the antisymmetry transform is 
{f q } = (f q {s) — e^ sxs ^(p q (s') . Such a specific feature originates from the generic 

permutation antisymmetry of the Fermi wave function of our many-electron system. We 
may therefore consider only 'antisymmetric' functions for which 



( Pq = {<Pq}/ 2 



(9) 



Our task is to find the energy of the eigenvector |X q ) and the 'wave function' (p q (s), assuming 
that the latter is regular and square integrable. If E q is the correlation part of the total CSFE 
energy (namely, Eqsfe = E v3iC +\gfiBB\ + hu c + E q ), then E q is found from— 



flED + flnT5^q \^)=E q \Xq). 



(10) 



Now we project this equation onto two basis states |p, q) and Qq— q 1 ) ? an d obtain two closed 
coupled equations 

(2i\g 1/2 
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and 

S oJ (q) + (2^)- 1 /^^( s) ( |Q 0Tq [H' 01 , Ql^Ql^] |0) = E q (12) 

s 

for E q and <p g (p). 

Next step is a routine treatment of Eqs. (11) and (12) in terms of calculation of commu- 
tators guided by commutation rules (3). In the q = case, which we immediately consider, 
the function v 2 o(p) depends only on the modulus of p. As a result we obtain— 



POO 

[E-S sw (q) -£ mp (<?)]y?o(<?) + / sds 

Jo 



K 1 (s,q)if (s) 



+ r^(l-co 8 [s x q])y?o(|q+s|) 

71 Jo 



(13) 



and 



1 I" 

E-£ HF = —J^ dpp 3 V(p)e-v 2 / 2 Vo (p) (14) 
we omit subscript in Eq), where 

„2 i poo 

g(q) = -±-V{q)e-« 2 ' 2 - — / p 3 V(p)e^ 2 / 2 J (pq)dp, (15) 



2y/2 2V2Jo 

K 1 (q,s) = ^e- s2 / 2 V(s)J (qs), and K 2 (s) = ^2-£) V{s)e- s2 l 2 (16) 

(0 in Eq. (13) is the angle between s and q). 

The problem has thus been integrable to yield in the thermodynamic limit a pair of 
coupled integral equations for one-dimensional function (fo{q) and the eigenvalue E. In 
order to solve this system we employ the method of expansion in orthogonal functions 



2N-1 



n=l,3,5,... 

These ip n = y/2L n (q 2 )e~ q2 ^ 2 with odd indexes of the Laguerre polynomials (f Q °° qdqip m ip n = 
5 mtn ) are chosen as a natural basis satisfying: (i) the property of integrability and ex- 
pected analytic and asymptotic features of <po(q); (ii) the antisymmetry condition (9). In 
other words, we change from the basis formed by the set of nonorthogonal double-exciton 
states |s, 0) = Qq^QqiJO) to a new set of basis states |DX, n) = {2N^)~ 1 ^ 2 ^ s i[) n (s)\s, 0) 
which are strictly orthogonal. Indeed, one can check by employing Eq. (3) and identity 
_i_ ^ s e^")' = / °° sds.J (qs)ip n (s) = -i/j n {q) that (m, DX|DX, n) = 8 m>n . The integer 
number N is dimensionality of this new double-exciton basis. 

After substitution of Eq. (17) into Eq. (14) the latter takes the form: E = F, where 

2JV-1 f00 

F = £ HF + — V A n dpp 3 V(p)e~ p2 / 2 Mp)- (18) 

V 2 n=l,3,5,... J ° 

Let us consider the ideal 2D case where V(q) — 1/q. (Here and below energy is mea- 
sured in units of e 2 / kIb-) After substitution of the expansion (18) into Eq. (13), further 
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multiplication by basis functions ip m (q) and integration ( f ...qdq) lead to the set of N linear 
algebraic equations with respect to A n . Finding A n for a given E and substituting them into 
Eq. (18), we obtain F(E). The condition F{E) = E yields the desired result E = E SF . 




0,0 0,2 0,4 0,6 0,8 

FIG. 1: Graphical solution of Eqs. (13) and (14). Intersection of the F = E straight line with the 
dotted line corresponds to the CSFE energy, -Esf~ 0.71 See text for details. 

Fig. 1 shows the result of calculations for N = 50. The lines which are restricted by 
vertical asymptotes reflect the result of calculation of F(E). Points of singularity E^\ at 
which F goes to infinity, are roots of the equation D N (E) = where D N is the determinant 
corresponding to the "left-side" of the set of equations for A n . By increasing iV we increase 
the order of equation D N (E) = 0, so that this has up to N real roots. Indeed, when 
observing the evolution of F{E) with increasing N, one finds that the number of singular 
points grows, and they become more densely placed. For N — > oo one could expect that a 
singular point appears within an arbitrarily small vicinity of every value E. Since all the 
vertical asymptotes E = are crossed by the straight line F = E (see Fig. 1), we come 
to the conclusion that for any E there is a singular solution of Eqs. (13) and (14). Such 
solutions with singular functions ip(q) form a band. The physical meaning of this result 
is quite transparent. Namely, the band corresponds to energy S sw (q) + £ mp (q) of unbound 
exciton pairs. Now we only consider the solution E = F(E), where the F = E line crosses a 
conventional envelope curve tracing the regions of regularity of ip determined by Eq. (17). 
Such regions at a finite N should be as distant as possible from the points of singularity, 
and we simply define them as the vicinities of "middle" points E® = \{E^ + E^ 1 ^). The 
envelope curve may obviously be defined as the line passing through the points [E , F(E )]. 
The intersection with the straight line F = E occurs at the only point stable with respect 
to evolution of this picture at N — > oo. This intersection point is readily seen in Fig. 1. 

Fig. 1 shows the build-up of singular points (vertical lines) with vanishing E and vice 
versa a certain rarefication of singularities in the vicinity of E$p. The former reflects growth 
of the density of states at the bottom of the exciton band whereas the latter is a usual effect 
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of the "levels' repulsion" . Note that the non-Hartree-Fock shift for the CSFE level is positive 
as compared to the value £'hf = 0.627. This is expected because the repulsion of the CSFE 
from the lower-lying crowded states of unbound excitons should be stronger than from the 
upper states having comparatively low density. At the same time, one can also see in Fig. 1 
some trend towards the concentration of singularity points at higher energies E. This 
is evidently a consequence of the density of states growth at the top of the exciton band. 

In general, the larger is N the more accurate is the calculation of <fio(q) and E, i.e. 
the envelope curve in Fig. 1 becomes discernible and may be drawn only at considerable 
N. At the same time the analysis reveals that the intersection point with the F = E 
line is rather stable and only weakly depends on N. This feature prompts us to consider 
the case N = 1 where double-exciton states mixed with Q Q j |0) are modelled by a single 
vector |DX, 1). Actually the N = 1 approximation for the problem determined by Eqs. 
(13), (14) and (17) is equivalent to a variational procedure for the trial double-mode state 
\X® M ) = QLJO) +AjjDX, 1^ where the correlation part of the excitation energy is found 
from equation 

t/vDMifr I yDM\ \ 
\ A l-"int|^-0 ) \ pint /-, q n 

W M |*o DM > ) ™ { ) 

{E^ c denotes the correlation part of the ground-state energy). After minor manipulations 

we find that this simple double-mode approximation (DMA) reduces our problem to the 

secular equation det\(E — £i)5ik + (1 — 5^)2)^1 =0 (indexes i and k are 1 or 2), where £i = 

J °°qdqV(q)e(q), £ 2 = £ HF , and D 12 = D 21 = J™ qdqV (q)d(q) with e = 2g 2 (l -q>)* e -W + 

|(g 2 -5g 2 +g 4 )e- < ? 2 -^(g 2 -4) 3 e- 3 ^/ 4 +(g 2 /2-2)e-" 2 / 2 and d = q 2 {q 2 -l)e~i 2 . Only the largest 

root of this secular equation has physical meaning. In the ideal 2D case we easily obtain the 

DMA correlation energy of the CSFE: E$f — 0.766. Comparing this result with Fig. 1 we 

conclude that even the DMA works rather well. 

Fig. 2 shows the CSFE correlation energy calculated within the DMA and employing 
the HF approximation, if the vertex V for a real 2DEG is defined as V = Ft(q)/q with the 
formfactor F b (q) = | (l + §) 3 8 + 9| + 3 (|) 2 Here b = b l B is a dimensionless parameter 
corresponding to dimensionless q. (bo is considered to be independent of the magnetic field.) 
It is seen that the non-HF shift of the CSFE energy, being about 15% in the strict 2D limit 
(i.e., in the b — ► oo case), becomes smaller (~ 5 — 6%) in real samples. This difference is 
not observable experimentally— Meanwhile, the DMA results are in good agreement with 
experimental data where the CSFE correlation energy is measured as a function of magnetic 
field, see inset in Fig. 2. The chosen value, bo = 0.213/nm, is quite consistent with the 
available wide quantum wells.— 

In conclusion, we note that preliminary analysis indicates that the non-HF shift should 
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FIG. 2: Main picture: DMA and HF shifts in dimensionless units against the form-factor parameter 
b. Inset: DMA shift against the magnetic field when b = 5.45 B~ x l 2 (b = 0.213 Ifl/nm, Ib in nm's, 
B in Teslas); symbols are experimental data for the 25 nm quantum wells.— 

be more substantial in the case of a fractional filling, e.g. at v— 1/3. Moreover, contrary 
to the single-mode approximation 8 shifting the energy to lower values as compared to the 
HF result, the approach taking into account the double-exciton component should lead to a 
considerable positive shift in the CSFE correlation energy. 
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1 If using ER, the relevant operators in our case are Hed= £ H a + Yl ^ab, where 

<z=0,0,l ab=00,01,01 

it n /v -1 / 2 nt 



Ha = \ £ q n?)>&(q) (QUQaaq- V S»aO, 

H ab = Y:^y{Q) [Mq)Mq)QUQMq+IMq)l 2 <W 6 (Ql^Qab^-N-^Ql 

and ^ritfH*^ 



L^ ria {q 2 /2)e~ q2 / 4 (Li\ is Laguerre polinomial). 



15 p or re f erence we wr ite out Eqs. (11) and (12) in the q 7^ case: 

[E q - £ sw (|q/2-p|) - £ mp (|q/2+p|)] ^(p)-^ {g q (p)} 
= (2ir)- 1 fds<p q (s) [[[/bo(|p-s|)-tT l(|q/2 + s 

+Uoi(|p - B|)e*(" x P)- - U 00 (\p - s|)e^ x ( s -P))-/ 2 



g i(pxs) 2 

U 01 (\p- 



s e 



i(qx(p-s))./2 



and 



E q = %(?) + ^ Jdpg*{p)<p q (p) 



with the 'free' term ^{<?g(p)}j where g q (p) 
and U nanb = V (q)h 2 ab {c\) (see notations of Ref. 



C/oi(|p + q/2|)e^ x ^/ 2 , Uoi=V(q)\h 01 (q)\ 2 



14). If q is chosen parallel to y, then <p q (p) is an 



even function with respect to the replacement p x — > —p x . The HF result £gy(<z) was calculated 
in Ref. 



